Deriving the cosmic ray spectrum from gamma-ray observations 
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A fundamental problem of cosmic ray (CR) physics is the determination of the average properties 
of Galactic CRs outside the Solar system. Starting from COS-B data in the 1980's, gamma-ray 
observations of molecular clouds in the Gould Belt above the Galactic plane have been used to 
deduce the Galactic CR energy spectrum. We reconsider this problem in view of the improved 
precision of observational data which in turn require a more precise treatment of photon production 
in proton-proton scatterings. We show that the spectral shape dN/dp oc p -2,85 of CR protons as 
determined by the PAMELA collaboration in the energy range 80 GeV < pc < 230 GeV is consistent 
with the photon spectra from molecular clouds observed with Fermi-LAT down to photon energies 
E ~ 1-2 GeV. Adding a break of the CR flux at 3 GeV, caused by a corresponding change of the 
diffusion coefficient, improves further the agreement in the energy range 0.2-3 GeV. 
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PACS numbers: 98.70.Sa, 13.60.Le, 98.38.Dq, 98.70.Rz 

I. INTRODUCTION 

The propagation of cosmic ray (CR) protons and nuclei 
with energies E/Z < 10 18 eV in the turbulent component 
of the Galactic magnetic field resembles a random walk 
and can be described in general by the diffusion approx- 
imation [l|. Therefore the Galactic disk should be filled 
with a well-mixed "sea" of CRs whose properties are sum- 
marized by the differential diffuse intensity 1(E) or the 
differential number density n(E) = An 1(E) /v. Exclud- 
ing the regions close to recent CR sources, the gradient 
V ln[rt(i?)] induced by the small current of CRs diffusing 
out of the disc and its extended CR halo is small. 

Most of our knowledge about Galactic CRs is obtained 
via observations in our local environment. Despite of the 
diffusive propagation of CRs, the CR intensity deduced 
locally may differ from the one averaged over the Galactic 
disc: At low energies, the influence of the Solar wind on 
measurements of the CR energy spectrum and the total 
CR energy density has to be corrected based on current 
understanding of the heliospheric modulation and direct 
CR measurements at different heliospheric distances and 
at different modulation levels. Clearly, such a correc- 
tion is model-dependent and can introduce uncertainties. 
Moreover, the Sun is close to a region with increased star- 
formation and thus supernova rate. Other reasons for 
local deviations include stochastic re-acceleration in the 
local interstellar turbulence or local sources as old super- 
nova shocks, winds and flares of massive stars. Therefore, 
the CR density close to the Solar system may deviate 
from the average Galactic one. 

A way to obtain independent information on the av- 
erage "sea" Galactic CRs is the observation of suitable 
molecular clouds far from CR accelerating regions 0: 
These clouds serve as a target for CRs producing gamma- 
rays mainly through decays of neutral pions created as 
secondaries in CR-gas collisions. Suitable clouds should 
be located away from the Galactic plane in order to test 
the "sea" CR spectrum, excluding the directions towards 



the inner and outer Galaxy. Assuming that gamma-ray 
production in hadronic interactions is sufficiently well- 
understood, the observed gamma-ray flux T^(E) from 
these clouds can be inverted to obtain the differential 
CR number density n(E). 

Previous works used observations of molecular clouds 
in the Gould Belt, in particular of Orion A and B, per- 
formed first by COS-B 0, then EGRET 0, Q, and 
most recently Fermi-LAT to derive the spectral shape 
of Galactic CRs. During this period, the quality of ex- 
perimental data has been hugely improved: On the ob- 
servational side, the data from Fermi-LAT have a much 
reduced error compared to its predecessor EGRET and 
extend now up to photon energies E 1 ~ 100 GeV, corre- 
sponding to typical energies of CR primaries E ~ 1 TeV. 
Moreover, the PAMELA collaboration determined the 
slope /3cr of the CR spectrum dN/dp oc p~^ CR with 
an accuracy of A/3cr = ±0.05 in the energy range 
80 GeV < pc < 230 GeV @. Thus the prediction of the 
secondary photon spectrum requires either similar precise 
photon fragmentation functions, or at least an estimate 
of their error. Finally, there are new results on photon 
yields from HERA Q as well as from LHC Q on the ac- 
celerator side, restricting theoretical models for photon 
fragmentation functions. 

In view of the improved precision of the experimental 
data, we reconsider this problem, paying special atten- 
tion to the treatment of photon production in proton- 
proton scatterings. We find that several commonly used 
paramctrisations for the photon fragmentation function 
as the ones of Refs. 0, [l(j deviate substantially from 
experimental data at high energies. These differences 
diminuish considering the photon yield produced by CRs 
with a power-law momentum distribution. In this case, 
wc find a relatively good agreement concerning the shape 
of the photon spectra, while the absolute photon yield 
differs by ~ 20%. As our main result, we show that the 
spectral shape dN/dp oc p~ 2 - 85 of CR protons as deter- 
mined by PAMELA in the energy range 80 GeV < pc < 
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230 GcV [6J is consistent with the photon spectra from 
molecular clouds observed by Fermi-LAT down to ener- 
gies E ~ 1-2 GeV. The agreement is further improved, 
if the CR spectrum exhibits a break around 3 GeV, as 
suggested by radio data [IT| . 

This work is structured as follows: We compare first in 
Sec. Inl several models used for the calculation of photon 
production in hadronic collisions to experimental data. 
We conclude that a combination of the parametrization of 
Ref. [H for nondiffractive processes below £" t hr = 50 GeV 
with the QGSJET-II model [l2j at higher energies gives 
a satisfactory description of experimental data. Then we 
calculate in Sec. IIIII the photon spectra expected from 
molecular clouds for a given CR flux. In the appendix, 
we describe the use of the photon and antiproton frag- 
mentation functions employed by us which are available 



from http : //sourcef orge .net/projects/ppf rag 





10 GeV 


100 GeV 


ITeV 


lOTeV 100 TeV 


QGSJET-II-04 


0.62 


0.64 


0.64 


0.70 0.79 


SIBYLL 2.1 




0.71 


0.84 


0.93 1.1 


Model of Ref. [9] 


0.57 


0.78 


0.73 


0.84 1.1 


ND model of [9] 


0.48 


0.53 


0.48 


0.49 0.69 



TABLE I: Predictions for the spectral-weighted moment 
Zry(Eo) (in mb) for photon production in pp collisions at dif- 
ferent laboratory energies Eq: Comparison of QGSJET-II-04, 
SIBYLL 2.1, the parametrization of Ref. 0, and the non- 
diffractive (ND) model of Ref. H; all for /3 C a = -2.85. 



In the next subsection, we compare the results derived 
using the "Kamae" parametrisation for the photon frag- 
mentation function given in Ref. @ to those obtained 
from the Monte Carlo simulations with the QGSJET-II- 
04 [HI and SIBYLL 2.1 [3 models, the latter having 



II. MODELS FOR PHOTON PRODUCTION 



been used as the basis for the parametrization of the 
photon fragmentation function of Ref. [l(| ■ 



High-energy photons can be produced both by CR pro- 
tons and electrons. In the latter case, inverse Comp- 
ton scattering on photons mainly from the cosmic mi- 
crowave background and bremsstrahlung are potentially 
contributing processes. In particular, bremsstrahlung 
was discussed as an important contribution to the to- 
tal observed gamma-ray spectrum from molecular clouds 
at energies below 100 MeV [f|. In this work, we restrict 
ourselves however to the energy range E 1 > 200 MeV ob- 
served by Fermi-LAT where bremsstrahlung can be ne- 
glected. For the density of molecular clouds, also inverse 
Compton scattering gives a negligible contribution rela- 
tive to photon production in CR-gas collisions. 

Photon production in hadronic collisions results 
mainly from decays of neutral pions produced as secon- 
daries. At sufficiently high energies, an additional though 
much smaller contribution comes from rj decays, while 
direct photon production is strongly suppressed and neg- 
ligible for astrophysical applications 1 . 

The photon yield in hadronic collisions can be calcu- 
lated using either numerical parametrisations or Monte 
Carlo simulations. The former are typically based on 
theoretically motivated or empirical scaling laws fitted 
to accelerator data. In general, they are well-suited for 
collisions at relatively low energies E < 50 GeV. In con- 
trast, Monte Carlo simulations are developed mainly for 
high-energy collisions and are based on a combination of 
non-perturbativc models and perturbative QCD. They 
treat both soft interactions and the hadronization of par- 
tons produced in (semi-) hard scattering processes. 



1 In contrast, the results of Ref. [l3| seem to indicate a significant 
contribution from direct photon production. This surprising re- 
sult is explained by the simple fact that the authors of Ref. [l3l | 
refer misleadingly to photons from r\ decays as direct photons. 



A. Comparison of the Kamae parametrisation and 
QGSJET-II 

We start by recalling some basic analytical formula 
before we discuss our numerical results for the various 
quantities characterizing photon production in hadronic 
collisions. Assuming a power-law cosmic ray spectrum, 
dNcR./dE = N E~^ CR , the resulting 7-ray flux may be 
written as 
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where xe = E^/E' and we have assumed -Emax ^ u 7 . 
Thus, any difference in the spectral shape between the 
parent CRs and the produced photons is introduced by 
the violation of Fcynman scaling, i.e. by the energy de- 
pendence of the spectral moment Z 7 . Such a dependence 
emerges because of i) the slow energy rise of the inelas- 
tic pp cross section ^^(E), ii) the relatively fast in- 
crease of the central rapidity plateau of secondary parti- 
cles 1/c™ g1 dijpp^ /dy\ y= o, and hi) the slow "softening" 
of the forward spectra of secondary mesons. While the 
first two effects result in an energy rise of Zj, the third 
one works in the opposite direction. 

For our purposes, it is more convenient to consider 
moments Z y (Eo) defined for a given energy Eq of the 
CR proton in the pp collision, 



Zj(E ) 



B ^jg^fg), (2) 



10 dx E 
In TablcO] we illustrate the energy dependence of Z 1 (Eq), 
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using /?cr = —2.85 as reported by the Pamela Col- 
laboration in the range 50 GeV < E < 200 GeV 
q|i comparing the predictions of the QGSJET-II-04 
Hi and SIBYLL 2.1 Q Monte Carlo generators 
to the results obtained using the parametrization of 
da pp ^'i{EQ,E 1 )/dE 1 from Ref. @. Clearly, the consid- 
ered models predict a quite different behavior of Zj(Eq) 
in the energy range of interest: While in the case 
of QGSJET-II the spectral moment is approximately 
energy-independent up to E Q ~ 1 TeV, Z^(Eq) has a 
relatively steep energy rise in the other two models. 

In order to decide which of these approaches provides 
a better description of photon production, we compare 
them next to data from accelerator experiments. We 
start by considering the respective results for photon 
spectra in proton-proton collisions at 205 GeV/c labora- 
tory momentum in Fig. [1] and for spectra of neutral pions 
and etas in pp collisions at p^ = 250 and 400 GeV/c in 
Figs. H and |3 In the case of QGSJET-II, the compari- 
son demonstrates a good overall agreement between the 
predictions and the experimental data. The Feynman 
Xf spectra of photons and neutral pions obtained with 
SIBYLL 2.1 are very similar to the QGSJET-II case at 
small Xf, but become substantially harder with increas- 
ing energy for xf > 0.2, see Figs. QJJD 

As the accuracy of the fixed-target data considered 
does not allow one to discriminate between these two 
trends, a valuable benchmark is provided by recent 
HERA measurements of photon production in the pro- 
ton fragmentation region for pr/* interaction at y^s = 319 
GeV Q: The observed forward energy spectra of gam- 
mas appear to be well-described by QGSJET-II, while 
being substantially softer than SIBYLL predictions. 

At first sight, this conclusion seems to be in a con- 
tradiction with the results of spectrometer studies of pp 
collisions at y/s = 7 TeV by the LHCf collaboration Q: 
The measured very forward photon spectra proved to 
be significantly harder than predicted by QGSJET-II. 
However, the preliminary results on the ir° production 
at y/s = 7 TeV by the same collaboration indicate that 
the latter discrepancy is likely due to a somewhat softer 
than observed transverse momentum distribution of neu- 

H 2 



tral pions in QGSJET-II 



The latter conjecture 



is supported also by the HERA data of Ref. Q which 
demonstrate that the pt-spectra of forward photons are 
better described by SIBYLL 2.1 than by QGSJET-II. 
We conclude that QGSJET-II agrees well with the ex- 
perimental data on the energy spectra of photons. The 
predicted pt distributions are somewhat too soft, but are 
anyway irrelevant for our problem. 

Though the considered Monte Carlo simulations are 
not designed to treat hadronic collisions for E < 50 



GeV, the QGSJET-II results appear to be relatively rea- 
sonable down to ~ 12GcV/c laboratory momentum, 
as demonstrated in Fig. [5j However, extrapolating the 
model to even smaller energies is meaningless, because 
the relevant physical processes, like resonance production 
and secondary Rcggcon exchanges, are not included. In- 
deed, as it can be clearly seen in Fig. El a t 8.8GeV/c lab- 
oratory momentum QGSJET-II tends to predict a harder 
photon spectrum than observed experimentally. 

In contrast, the situation with the model of Ref. Q 
appears to be quite different: This parametrisation de- 
scribes the experimental data quite well up to Eq ~ 
10 GeV, cf. Figs. [5] and U Figures □ and H show that 
the photon spectra predicted by the Kamae parametriza- 
tion become at higher energies much harder than those of 
QGSJET-II which, as we have argued above, agree with 
HERA data. This explains the larger values for Zy(Eo) 
obtained in that case in the energy range Eq > 100 GeV, 
see Table 1. The discrepancy between the parametriza- 
tion of Ref. [j| and the data has its origin in the somewhat 
oversimplified treatment of diffractive particle produc- 
tion in the underlying model [Tsj ] which utilises a cluster- 
like hadronization procedure: The energy of the diffrac- 
tive state is distributed more or less uniformly between 
the pions produced, neglecting the leading baryon effect. 
In reality, a large part of high energy diffraction corre- 
sponds to the creation of high mass states which are de- 
scribed by the Pomeron contribution and correspond to 
multiperipheral kinematics of particle production 3 [l7| . 
As a consequence, forward spectra of neutral pions and 
etas in high mass diffractive processes resemble the ones 
of the usual nondiffractive collisions. 

The complementarity of the two models motivated 
us to combine the QGSJET-II description with the 
parametrisation of Ref . Q: While we use QGSJET-II for 
photon production at relatively high energies Eq > E t hr, 
with 10 < Ethr < 50 GeV, at lower energies we apply 
the parametrization of Ref. [i[ restricted to nondiffrac- 
tive processes. 



III. GAMMA-RAY SPECTRA 

We explore now the consequences of our new photon 
fragmentation function. Let us remind first that diffu- 
sive shock acceleration predicts a power-law in momen- 
tum, dN/dp cx p _/3cR , cf. Ref. [l]|, while the interstellar 
propagation is also rigidity-dependent. It is, therefore, 
natural to use a power-law in rigidity to fit CR data. 
Indeed, PAMELA p and He spectra and their ratio can 
be well-described with a power-law in rigidity down to 
~ 20 GV, below which the interstellar spectrum is signif- 
icantly modified by the heliospheric modulation. 



2 Fixed angle spectrometer measurements of forward particle spec- 
tra are very sensitive both to the respective Feynman x and pt 
distributions, the two variables being related to each other in the 
cm. frame as pt = i? Q bs */s/2. 



3 Even more complicated diffractive final states are produced at 
much higher energies Il2t 1 1SH . 
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FIG. 1: Photon production cross sections in the center of mass (cm.) frame for proton-proton collisions at 205 GeV/c: 
Feynman x spectrum (left), rapidity distribution (middle), and transverse momentum distribution (right). Calculations with 
QGSJET-II-04 (red solid), SIBYLL 2.1 (blue dot-dashed), parametrization of Ref. (green dashed), and ND part (black 
dotted line) of the latter are shown together with experimental data from Ref. [2lj . 




FIG. 2: Feynman x spectrum (left) and transverse momentum distribution (right) of neutral pions in the cm. frame for pp 
collisions at 250 GeV/c as calculated using QGSJET-II-04 (red solid) and SIBYLL 2.1 (blue dot-dashed) compared to the data 
from Ref. [H. 



In Fig. we compare the gamma-ray spectrum pro- 
duced by cosmic rays with dN/dp oc p~ 2 using QGSJET- 
II, the parametrization of Ref. or using only the 
nondiffractive part of the latter. The photon spectrum 
obtained using the Kamac parametrization rises much 
quicker with energy in the 1-10 GeV range than the one 
obtained with QGSJET-II, as a result how the diffrac- 
tion is modeled in Ref. [l(| , while at higher energies the 
obtained spectra have similar shapes, with ~ 20% dif- 
ference in the normalization. On the other hand, the 
results obtained using the nondiffractive part of the Ka- 
mae parametrization agree well with QGSJET-II up to 



few tens of GeV - compatible with the similarity of the 
photon spectra for pp collisons in the two models in the 
energy range E p ~ 10 — 200 GeV seen 4 in Figs. Q] and |U 
Therefore, we find it natural to match in the following 
the results of the QGSJET-II model for photon spectra 
to the parametrization of Ref. Q for nondiffractive pro- 
cesses at -Ethr = 50 GeV. 



4 The difference at 1 — 3 GeV is caused by the extrapolation of the 
QGSJET-II model outside its working range: Below 10 GeV lab. 
energy, it predicts too hard photon spectra, as shown in Fig. \E\ 
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FIG. 3: Feynman x spectra for ty° (left) and r\ (right) production in the cm. frame for pp collisions at 400 GeV/c as calculated 
using QGSJET-II-04 (solid) and SIBYLL 2.1 (dot-dashed) compared to the data from Ref. [H. 




FIG. 4: Calculated energy distributions of photons in the laboratory frame for pp collisions at 50 GeV, 5 TeV, and 500 TeV; 
abbreviations for the lines are the same as in Fig. [T] 



As a further application of our new photon fragmen- 
tation functions, we consider the problem how the cos- 
mic ray spectrum can be derived from gamma-ray ob- 
servations. The Fermi LAT collaboration has studied 
this question in great detail and using a variety of meth- 
ods [28l - l3l| : For instance, the diffuse Galactic gamma- 
ray emission was compared in Ref. [28| to GALPROP 
models and an overall agreement of ~ 15% of data and 
models was found. In Ref. [2{|, the analysis of the 
diffuse gamma-ray emissivity was constrained to well- 
defined segments of the Local and the Perseus arms, 
deriving thereby constraints on the cosmic ray density 
gradient. Using observations in the mid-latitude region 
in the third quadrant, the Fermi-LAT callaboration con- 
cluded in Ref. [3(| that the CR spectrum derived agrees 



with the locally measured one within 10%. 

More recently, the authors of Ref. [32j used Fermi-LAT 
observations of nearby molecular clouds to deduce the 
energy spectrum of Galactic sea CRs. The photon flux 
deduced in Ref. [32[ is shown in Fig. [5] as red error-bars 
together with the photon spectrum (dashed green line) 
derived by us assuming a cosmic ray spectrum charac- 
terized by the slope /3cr = —2.85, and using the combi- 
nation of QGS JET-II results for E p > 50 GeV with the 
parametrization of Ref. @ for nondiffractive processes 
at lower energies. The predicted photon flux agrees well 
with the data at energies E y > 1 GeV. The remaining 
discrepancy at lower energies is at most at 30% level and 
may be partly related to uncertainties in the description 
of photon production at Eq < 50 GeV 0] . A more im- 
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FIG. 5: Longitudinal momentum spectrum (left), rapidity distribution (middle), and transverse momentum distribution (right) 
of photons in the cm. frame for pp collisions at 12.4 GeV/c as calculated using QGSJET-II-04 (solid lines) compared to 
experimental data from Ref. [24| . Dashed and dot-dashed lines - parametrizations of Ref. [t| for photon production in inelastic 
and nondiffractive pp collisions, respectively. 




FIG. 6: Calculated Feynman x (left) and rapidity (right) spectra of photons in the cm. frame for pp collisions at 8.8 GeV/c 
compared to experimental data from Ref. [25| ; abbreviations for the lines are the same as in Fig. [5] 



portant reason for the difference at energies E 1 < 1 GeV 
may be the energy-dependence of the diffusion coeffi- 
cient, which is expected to change at 2-3 GeV [IJ H3] : 
For instance, Ref. [ll| found fitting Galactic synchrotron 
data that the diffusion coefficient reaches a minimum 
at 2-3 GeV, increasing as p^ 1 or even faster at lower 
momenta p. To illustrate the latter point, we calcu- 
lated the diffuse gamma-ray flux for the case of a broken 
power-law CR spectrum: with the slope /?cr = —2.85 for 
Pgr > 3GeV/c and /3 C r = -1 for Pc r < 3GeV/c. The 
result, plotted in Fig. [8] as the solid blue line, matches 
well the observations in the whole energy range consid- 
ered. 



Next we examine, if a spectral break at higher en- 
ergies is also consistent with the Fermi-LAT data. In 
particular, a broken power-law with high-energy slope 
/?cr — —3, low-energy slope /3cr — —1-9, and break- 
energy at T = E p — m = 9 GeV, or E p ~ 10 GeV, has 
been obtained in Ref. [32[ as best-fit to the Fermi-LAT 
data. In Fig. |9l we show as blue solid line the gamma- 
ray flux corresponding to a CR spectrum with the slope 
/3 C R = -3 for pcr > fOGeV/c and f3 CR = -2.4 for 
Pcr < lOGeV/c which agrees clearly with the data also 
well. Note that the sharper spectral break obtained in 
Ref. [H], with the slope /3 C r = -1.9 for T < 9 GeV, 
is caused by the use of the power-law CR spectra with 
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FIG. 7: Gamma-ray flux E 2 T{E) produced by cosmic ray 
protons with a spectrum dN/dp cc p~ 2 as calculated using 
QGSJET-II (solid red line), the Kamae parametrization [S| 
(dashed green line), or the nondiffractive part of the latter 
(dotted black line). 
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FIG. 8: Gamma-ray flux E 2 J r {E) as calculated using the 
combination of QGSJET-II and the ND model of Ref. |<J 
(£thr = 50 GeV) for a single power-law CR flux with the slope 
/3cr = —2.85 (dashed green line) and for a broken power-law 
with /3cr = —2.85 for pgr > 3 GeV/c and /?cr = —1 for 
Pcr < 3 GeV/c (solid blue line). The gamma-ray spectrum 
derived in Ref. [32| from Fermi-LAT data is shown as points 
with error-bars. 
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FIG. 9: Gamma-ray flux E 2 J r (E) for a broken power-law with 
Pen = —3 for pen. > lOGeV/c and ^ch = —2.4 for pen < 
lOGeV/c as calculated using the combination of QGSJET-II 
and the ND model of Ref. Q for E thr = 50 GeV (blue solid 
line) and based on the Kamae parametrization alone (green 
dashed line); All fluxes are normalized to coincide at 1 GeV. 
The gamma-ray spectrum derived in Ref. [S^l from Fermi- 
LAT data is shown as points with error-bars. 



the spectral shape. 

Although the CR spectrum derived in Ref. (HJ is con- 
sistent with the data, we consider it as less attractive: 
First, this solution requires additionally to the break at 
9 GeV another one around 200 GeV, where the transition 
to the slope measured by PAMELA should take place. 
Second, both break energies have no obvious physical 
reason. In particular, it is surprising that the CR spec- 
trum measured inside the Solar system differs substan- 
tially from the one of Galactic "sea" CRs up to 200 GeV. 
In contrast, the spectrum shown earlier in Fig. |8]has a 
single break at an energy where a change of the diffu- 
sion coefficient is expected [Tl], [2(| and a single expo- 
nent /?cr = —2.85 in agreement with PAMELA mea- 
surements. 



IV. CONCLUSIONS 



respect to the kinetic energy in Ref. [32J, leading to a 
substantial enhancement of the region of small pgr- To 
test the sensitivity of these results to the photon frag- 
mentation function used, we repeat the calculation using 
the Kamae parametrization: We obtained a very similar 
photon spectral shape, as illustrated by the dashed green 
line in Fig. O though with a 20% higher flux. For an eas- 
ier comparison of the shape, we normalize all the fluxes 
to coincide at 1 GeV. As we have shown above, the differ- 
ence between the various fragmentation model manifest 
themselves mainly in the absolute photon yield, not in 



We have reconsidered the problem of determining the 
average properties of Galactic CRs using gamma-ray ob- 
servations of molecular clouds. The largely improved 
quality of the observational data requires a careful treat- 
ment of the photon fragmentation function. Comparing 
photon fragmentation functions calculated in different 
approaches at high energies, we have argued that a com- 
bination of the Kamae parametrisation and QGSJET-II 
provides the most reliable results. As our main result, we 
obtained that the spectral shape of CR protons as deter- 
mined by PAMELA is consistent with the photon spectra 
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from molecular clouds observed by Fermi-LAT down to 
energies E ~ 1-2 GeV. The agreement is improved fur- 
ther, if the CR spectrum exhibits a break around 3 GeV. 
This gives additional evidence for a change of the diffu- 
sion coefficient around 3 GeV, which was previously sug- 
gested on theoretical grounds [2(| and supported by ob- 
servations fllj . 



E s particles the inclusive spectrum E s da(E p , E s )/dE s 
is then obtained via an interpolation between the tab- 
ulated values. In the case of pp collisions, the in- 
terpolation between the non-diffractive part of Ka- 
mae's parametrisation and QGSJET-II used in this 
work is provided. All procedures are available at 
http : / / sourcef orge .net /pro j ects/ppf rag 
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Appendix A: Using the fragmentation functions 

The photon fragmentation functions for proton-proton, 
proton-helium, and helium-proton collisions have been 
tabulated based on QGSJET-II simulations for a set of 
incident proton energies E p between 10 and 10 s GeV. 
For given energies of the primary E p and secondary 



Appendix B: Comparison of the antiproton 
fragmentation function with data 

In addition to photon production, we performed a com- 
parison of model predictions for antiproton spectra in 
proton-proton collisions with available data, being mo- 
tivated by the importance of those spectra in particular 
for dark matter searches. 

As one can see from Figs. [TU] and [TT1 prc-LHC data 
are somewhat uncertain concerning the antiproton yield 
in pp collisions, which partly explains the vast differences 
between various model predictions. 

An important benchmark on the energy-dependence of 
the antiproton production is provided by the recent data 
obtained at LHC [H|, as illustrated in Fig. [T2] The an- 
tiproton yield in QGSJET-II-04 has been adjusted using 
these data. 

Fragmentation functions for the production of an- 
tiprotons (plus antineutrons) in proton-proton, proton- 
helium, and helium-proton collisions arc also available at 
|http : / / sourcef orge .net /pro j ects/ppf rag. 
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FIG. 12: Transverse momentrum spectra of antiprotons at 
central rapidities (\y\ < 0.5) in the cm. frame at y^s = 900 
GeV (open squares) and 7 TeV (full squares) as measured by 
the ALICE Collaboration compared to calculations with 
QGSJET-II-04 (red solid) and SIBYLL 2.1 (blue dot-dashed). 



